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Abstract. There are a number of approaches to testing General Relativity (GR) on linear scales 

using parameterized frameworks for modifying cosmological perturbation theory. It is sometimes 
assumed that the details of any given parameterization are unimportant if one uses it as a diagnostic 
for deviations from GR. In this brief report we argue that this is not necessarily so. First we show 
that adopting alternative combinations of modifications to the field equations significantly changes 
the constraints that one obtains. In addition, we show that using a parameterization with insufficient 
freedom significantly tightens the apparent theoretical constraints. Fundamentally we argue that it 
is almost never appropriate to consider modifications to the perturbed Einstein equations as being 
constraints on the effective gravitational constant, for example, in the same sense that solar system 
constraints are. The only consistent modifications are either those that grant near-total freedom, as 
in decomposition methods, or ones which map directly to a particular part of theory space. 
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1 Introduction 

The availabihty of high-precision cosmological data has made it possible to test General Relativity 
(GR) on cosmological scales. In particular, in the absence of a conclusive theoretical explanation for 
the dark energy phenomenon, modifications and extensions to that venerable theory remain a viable 
possibility. 

Since the process of comparing a new theory of gravity to cosmological data can be complicated, 
any test statistic that can be used to quickly reject a theory is a useful one. In this spirit, inspired by 
the Parameterized Post- Newtonian (PPN) formalism [1-4], a similar approach has been developed, 
which describes the deviations of a theory from the standard evolution of a perturbed GR universe 
[5-22]. Unlike the PPN case, however, such a parameterization is not unique; on cosmological scales 
there is freedom to modify evolution over the whole of cosmic time and over a range of scales. One 
is therefore led to constrain functions of space and time rather than just a handful of numbers in the 
weak-field metric. There are then two levels of structure that must be imposed on the parameterized 
field equations: 

Level 1: One must choose the way in which free functions are used to extend the field equations. 
For example, one could use free functions to rescale terms that exist in GR, or allow new terms 
to be present. 

Level 2: When implementing the parameterized field equations in numerical codes one is forced to 
choose a sensible ansatz for the time- and scale-dependence of the free functions. Typically 
these ansatzes are motivated by the observation that the dark energy-like sector only dominates 
the energy density of the universe at late cosmological times. See [23] for an enumeration of 
such ansatzes. 

The choice usually made at level 1 is to introduce two free functions directly into the modified Poisson 
equation and the 'slip relation' (the transverse, traceless component of the field equations), as these 
are the expressions relevant to observables such as weak lensing of galaxies and measures of structure 
growth [16, 17] . One of these free functions acts as an effective rescaling of Newton's gravitational 
constant, whilst the other is defined as the ratio of the two potentials that describe the perturbed 
metric in the Conformal Newtonian gauge. We will refer to this approach as 'phenomenologically- 
based'. It is not the case that particular choices of the two free functions are designed to reproduce the 
field equations of specific modified gravity theories (except in a handful of special cases). Rather, in 
the phenomenological approach, the free functions should be considered as indicator flags for non-GR 
behaviour. 
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A disadvantage of phenomenological-type parameterizations is their tendency to obscure which 
regions of theory space they correspond to, because there is no direct mapping between the free 
functions and the parameters of a specific model. It is therefore difficuh to translate constraints on 
the two free functions into constraints on a given theory of modified gravity. 

We have advocated an alternative approach: directly specifying quantities needed to describe a 
4x4 tensor of scalar modifications to the linearly perturbed Einstein equations [18-20]. We can then 
derive the corresponding Poisson equation and slip relation, which will contain components of this 
new tensor. This has been dubbed the 'Parameterized Post-Friedmann' approach (PPF), and it leads 
to a level- 1 structure (as classified above) which is different to that of the phenomenological approach. 

Arguments can be made in favour of both the phenomenologically-based and PPF approaches. 
For large classes of theories the two strategies become equivalent at intermediate distance scales where 
the time derivatives of perturbations can be neglected in comparison with their spatial derivatives. 
However, this approximation is not applicable on scales comparable to the horizon size. 

Whilst such large scales cannot be probed directly with galaxy surveys, they are nonetheless 
important for accurate calculation of the Integrated Sachs- Wolfe (ISW) effect and large-scale matter 
power spectrum. These quantities are frequently computed using Einstein-Boltzmann solvers such as 
Camb [24] or Class [25] which evolve perturbations through horizon crossing. 

The purpose of this brief report is to illustrate that the choice made for the level- 1 structure of 
a parameterization has a significant influence on the constraints obtained on deviations from GR - a 
caveat that is commonly forgotten. 

We demonstrate the differences between two possible parameterizations by generating some of 
the simpler perturbative observables - the CMB power spectra, in particular the ISW effect, matter 
power spectra, and the growth function f{z). In §4.2 we further show that the choice made for the 
level-2 structure has an equally important influence on the constraints obtained. 

2 Theory 

In this section we describe the two parameterizations that we apply in this paper. Our goal here is 
not to advocate for either of them, but rather to highlight the fact that they lead to different results. 
Parameterization A described below is derived as a special case of the general PPF form detailed in 
[26] . However, for the purposes of this paper we can ignore its origin and simply regard the resulting 
field equations as another example of a phenomenological-type parameterization - one with different 
choices made for the level-1 structure. 

2.1 Parameterization A 

In [18] a parameterization was proposed by writing the modifications as an additional tensor to the 
Einstein equations of GR: 

G^, = SnGoa^Tl^t + a'U^, (2.1) 

where the stress-energy tensor T^l contains all the standard cosmologically-relevant fluids and the 
tensor Uf^i, may contain metric, matter and additional field degrees of freedom coming from a theory 
of gravity. In writing (2.1) it is assumed that all known matter fields which are part of T^l couple 
to the same metric g^^, and that G^^ is the Einstein tensor of that same metric. This ensures that 
the stress-energy tensor of matter obeys its usual conservation equations, and hence U^^ is separately 
conserved. 

The formalism proceeds by parameterizing around the linearly perturbed version of equation 
(2.1). To enable a direct comparison with parameterization B below we will specialise to the case of 
purely metric theories, that is, those for which the action is constructed from curvature invariants 
(e.g. f[K) gravity) or non-local invariants (e.g. [27]). 
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The most general perturbations of [/^^ that are allowed in a second-order metric-only theory are 
as follows [18-20, 26]: 

-a^5U° = kBo^ 
a^SUi = k'^Co^ + fcCid 

a^SU'j ^ Dn4>+^^ (2.2) 

where dots denote derivatives with respect to conformal time rj and k is the Fourier wavenumber. 
The gauge- invariant metric perturbation $ reduces to the curvature perturbation $ in the Conformal 
Newtonian (CN) gauge. In our conventions the CN gauge is defined by: 

ds^ = -a^{l + 2^)dTf + a^{l - 2$)df2 (2.3) 

Further below we introduce a second gauge-invariant metric perturbation, ^ which reduces to ^' in 
the CN gauge. 

The coefficients Aq . . . Di appearing above are functions of background quantities, Aq — Ao{k, rj) 
etc., and the factors of k are chosen such that Aq . . . Di arc dimensionless. However, these functions 
are not all independent. Perturbations of the Bianchi identity [//f.^ = yield a set of additional 
constraint equations which can be used to reduce the six free functions in equations (2.2) down to 
just two. One of these free functions is defined to be: 

The modified Poisson equation can then be written (in Fourier space): 

- fc2$ = 47r-^a2pA = inGasa^pA (2.5) 

with Go denoting the canonical value of Newton's gravitational constant as measured by a Cavendish 
experiment on the Earth, and where the gauge-invariant comoving density perturbation pA is a 
summation over all conventional cosmologically-relevant fluids. The combination Gcfr = Gq (1 — .9)"^ 
plays the role of a modified Newton's constant. 

We choose the second free function to be Dq, which we will hereafter relabel as C = C{k,Ti) to 
distinguish it from its appearance in the more general format of equations (2.2). ( appears in the 
'slip' relation between the potentials <t and (which are equal to one another in GR): 

$-* = 87rGoa2(p + P)S + C3' + :|:$ (2.6) 

The anisotropic stress of conventional matter S is negligible after the radiation-dominated era. 

The key point here is that in parameterization A the slip relation and the modified Newton's 
constant are not independent, as the function g appears in both. This special case of the PPF 
formalism does not capture the behaviour of many modified gravity theories, because we have not 
allowed for additional degrees of freedom to appear [20, 26]. It is, however, directly comparable to a 
phenomenological-type parameterization. Hence we will regard eqns.(2.5) and (2.6) simply as possible 
alternatives to eqns.(2.9) and (2.10). 

2.2 Parameterization B 

If one steps back and looks at the structure of the evolution equations for cosmological perturba- 
tions, one finds that it is enough to work with only two of the four Einstein field equations: the 
Newton-Poisson equation and the slip equation. From a practical point of view, the only observ- 
able modifications to gravity (at least at the perturbative level) will be modifications to the these 
equations, which we can write in the following form 

- fc2$ - inGoo^pA = Fi (2.7) 
(* - $) + 87rGoa2(p + P)I] = F2 (2.8) 
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where Fi and F2 are arbitrary functions of time and space. A simple ansatz (from the dimensional 
point of view) is that Fi = ak^^ and F2 = — Reorganizing the equation we find that they can be 
rewritten as 



Parameterization B is defined in the CN gauge, so that $ and replace $ and in eqns.(2.5) and 
(2.6). This formulation is equivalent up to small corrections to that used in [16]; their anisotropic 
stress is related to ours by k^'S — a. This parameterization suggests that, once the anisotropic stress 
has become negligible, a theory of modified gravity could potentially modify one of equations (2.9) or 
(2.10) whilst leaving the other unchanged. This behaviour does not arise analytically from theories 
below third order [20]; but given enough freedom in the functional ansatz the numerical behaviour 
that occurs in any theory can be realised with this parameterization. (We note that a handful 
of parameterizations that are instead optimized for weak lensing analysis have also been suggested 



To demonstrate the differences that the alternative parameterization types can generate, we will 
use the same ansatz for the free functions Gcff and ^ in both parameterizations A and B. Since we are 
focusing on modifications to gravity associated with dark energy, an appropriate parameterization is 
an expansion in JIa- We therefore use a polynomial expansion up to third order in Q\. 

3 Spectra 

We have modified a version of the Boltzmann code Camb to use the new equations specified in section 
2. We use variants which include both parameterizations A and B. In this section we show and analyze 
CMB and LSS spectra generated by these models. We will show that the two variants, even when set 
up with equivalently parameterized functions, generate considerably different spectra. 

Figure 1 shows CMB spectra from parameterization A with variations to the linear term. Be- 
cause we have modified the late-time behaviour of gravity whilst leaving the era before dark energy 
domination unchanged, the small-scale CMB (which depends mainly on what happens at recombina- 
tion) is unaffected, but the large-scale features caused by the ISW effect vary. It can be seen that 
strengthening gravity by a small amount decreases the power spectrum of the ISW effect by counter- 
acting the suppressive effect of A on structure formation, whereas weakening GcS increases the ISW 
effect. This difference is significant because the ISW effect is only weakly present observationally. 

The behaviour of the spectra for a simple linear variation in Gcs shows the same general trends for 
model B as those for model A shown in figure 1. The most important difference between the theories, 
however, appears when the two free functions interact to allow cancellations. This is discussed in 
detail below. 

Examining the plot of linear variations in Gcff alone, one might be tempted to conclude that 
the perturbative strength of gravity, as measured by the Geff parameter, is rather well constrained 
using just CMB data. We wish to emphasise in this paper that this is not so. First, the various ways 
one can construct a theory mean that direct interpretation of the parameter, except on small scales, 
is not unique. Secondly and probably more importantly, the linear variation alone masks a richer 
phenomenology, which we will discuss more in the next section. 

Figure 2 shows the growth function, f{z) = scaled by erg, which is measured by, for 

example, redshift-space distortion experiments. Results from the recent WiggleZ survey [28] and 
other data [29-32] are also shown. These results provide weaker constraints than the ISW, but near- 
future experiments should reach sensitivity levels good enough for precision tests of late-time gravity. 

We do not display the matter power spectra because they are very mildly affected by the modi- 
fications of gravity, except on very large scales far beyond the reach of galaxy surveys. 

Figures 1 and 2 demonstrate that relatively small modifications of about 20% to the effective 
Newton's constant in the Poisson equation can have effects that are easily detectable at large scales in 
the CMB. But it is the combined effects of Geff and ^ that are of most interest. Figure 3 demonstrates 



fc2$ ^ 47rGofia^pA 

* = -87rGoa2(p + P)S + (1 



0* 



(2.9) 
(2.10) 



[11, 13, 15]). 
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how this combination can cancel even quite extreme individual effects. The curves in this plot have 
Goff/Go = f + 5.3r2A - 1-1 f^A + l-3f^A' and C = 4.8f7A - f.Of^^ + 6.61)^. The other cosmological 
parameters are within normal ranges, though not identical to those in figure f . There is a near-total 
cancellation of the ISW effect in model A, but model B with the same parameters is completely ruled 
out. There do exist alternative parameter combinations in model B with the same cancellation, but 
the fact that they are very different from those in model A is what we want to highlight here. 

We are not, of course, advocating these models as actual cosmologies - they simply demonstrate 
that a vast region of parameter space can be consistent with CMB data. Since other upcoming data 
sets like lensing and redshift-space distortions also constrain (different) combinations of the potentials 
there should be equivalent degeneracies present there. 

The importance of other data sets in breaking the degeneracy in the CMB is illustrated in figure 
4. That plot shows the growth rates (Jgf(z) for the extreme models described above. The behavior 
of the extreme case is once again very different between the two parameterizations. 

4 Constraints 

4.1 Parameter Estimation 

Now that we have demonstrated that the different parameterization approaches can yield considerably 
different results for power spectra, we can check how these differences go forward into parameter 
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Figure 3. A demonstration of how extreme changes to the effective gravitational constant can be permitted 
by the CMB when they can be counteracted by a significant gravitational slip (see text for form of extreme 
curves). The cancellation is absent for these parameters in model B. This delicate difference underscores the 
need for a rigorous treatment of modified functions. 
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Figure 4. The growth rates /(z) corresponding to the CMB spectra shown in figure 3. 



estimation. Note that at this stage we are not attempting to find the tightest possible constraints 
on these parameterizations using all available data. In particular we are not using data from weak 
lensing, galaxy-ISW correlations, or growth rates which can provide the strongest constraints. 

We are rather intending to demonstrate that the different schemes applied to the same data can 
produce significantly different constraints, even on quantities we usually regard as rather physical - 
such as the the effective Newton's constant controlling the growth of structures under gravity. We 
hope to show that headline constraints on the two free functions of parameterized approaches should 
be taken with a degree of caution. 

The plots in this section are generated from Monte-Carlo Markov Chains running parameteriza- 
tions A (equations (2.5) and (2.6)) and B (equations (2.9) and (2.10)). In each case, as in the previous 
section, the same ansatz is applied to the free functions Gcs(z) and C,(^z). The same constraining data 
is used in each case: the 7- year WMAP CMB data [33], the SDSS DR7 matter power spectrum [34], a 
prior i/o = 73.8 ± 2.4 [35], the BBN constraint l^fc/i^ ^ 0.022 ± 0.002 [36], and the Union2 Supernova 
la data [37]. ^ . With these data sets the strongest constraining power comes from the ISW effect (or 

^It is not generally correct to use published supernovae constraints directly when using modified gravitational physics, 
since the calibration factors applied to them are cosmology-dependent; it is only because we leave our background 



- 6 - 



rather lack thereof) in the large-scale CMB temperature power spectrum. 
4.2 InsufRcient ansatzes 

The simplest ambiguity we can consider when constraining the free functions can arise if we use what 
might be termed an incomplete parameterization - one where the free functions are not free enough to 
explore the available parameter space. One example of this would be our restriction here to functions 
only of z, with no scale-dependence. 

Figure 5 illustrates another example of an incomplete parameterization. In that figure we show 
how the constraints on the theory depend on the number of terms in the polynomial expansions of 
Gcff(^^A) and C(riA)- The larger contours use cubic models for the evolution, and the much smaller 
ones use only a linear term in the expansion. Both curves are for parameterization A, and the contours 
are for 68% and 95% probability mass. 

Of course we do not claim that the cubic model is general enough to describe the free functions 
sufficiently - it is the comparison that we highlight. Whenever contours are presented for the free 
functions of some parameterized model one should bear in mind that adding only a little freedom can 
change the areas of the error ellipses by a large factor. Hence all constraints of this form should be 
taken with a pinch of salt. 




Figure 5. Constraints on the free functions, evaluated at redshift zero, in a model restricted to linear evolution 
in (filled contours) and one with cubic evolution (line contours). The apparent constraints on the two 
parameters, for example, are misleadingiy small in the more restrictive parameterization. 



4.3 Variant parameterization 

The second result of this section is shown in figure 6; it demonstrates the parameterization-dependence 
of the constraints, independent of the ansatz. Most clear is that the free functions are significantly 
less correlated in model A. In the model B there is a simple and straightforward impact on the term 
whose derivative is the ISW source: 

(4.1) 



$ + *(X (C-2)Geff-^pA 
which yields a multiplicative correlation, whereas case A does not have this simple relation: 



$ + * OC (Geff - Go) 

+ (C - 3)Gcff + 



pA 



(4.2) 



1 - 



Go 
Gcff 



Gcff + Go 



fc2 



pA 



evolution unchanged from GR that is is possible here. 
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The edge of the curve for parameterization A in figure 6 cuts off rather sharply just below 
Gcff — 1. This is due to the impact of Gcff and C on the growth of matter density perturbations, and 
the restrictions that they should obey in order to reproduce the observed amount of structure in the 
universe today. During the matter-dominated epoch density perturbations on subhorizon scales grow 
as (5 cx a^' , where in parameterization A [38] : 

In parameterization B the exponent is different: 

PB = -l + \f^24^i^) (4.4) 

One can see that in the GR limit (C = 0, Gcff = Gq) we recover pa = Pb = 1, in agreement with 
the standard result for an Einstein-de Sitter universe. Equation (4.3) implies that in a universe with 
C = matter perturbations fail to grow during the matter-dominated epoch for values of Gcs/Gq < 
0.5. Hence we should expect this region of parameter space to be disfavoured by the matter power 
spectrum. The analogous boundaries in parameterization B will be different, which may explain 
the different shapes of the fall-offs of the distributions in figure 7. However, we ought to remember 
that equations (4.3) and (4.4) apply only in a pure Einstcin-de Sitter setting. The boundaries in 
parameter space that they imply may not be obeyed rigidly in the real universe, due to growth during 
the radiation- and A-dominated eras. 

For small deviations from GR Pa ~ Pb- This is because in EdS universes $ = so the slip 
equations (2.6) and (2.10) are the same. This does not apply when A becomes important. 

It is in the tails of the distributions where the difference between the two parameterizations can 
be most stark - this will be particularly important when we are trying to decide if GR is threatened 
by evidence of deviations. The histograms in figure 7 demonstrate this most clearly. In addition, the 
spectral index is much less constrained in model A than in either GR or model B, because in that 
model we can more freely cary Gcff at low redshift; so the matter power spectrum amplitude is not 
as constraining. 

3 
2 
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G„a(z = 0)/Go 

Figure 6. Joint constraints on the slip parameter <^ and Gcff at z = 0, for parameterization A (black lines) 
and B (filled green). In both cases 68% and 95% contours are shown. 



5 Discussion 

The first, and uncontroversial, issue we have noted here is the significant effect that using insufficiently 
free functions has on the constraints one obtains on modifications to the slip and Poisson relations. 




- 8 - 



2 4 6 8 -3-2-1 12 3 
Gefi(2 = 0)/Go = 0) 



Figure 7. Likelihoods of the effective gravitational and slip parameters at z = 0, for parameterization A 
(thick blue) and B (thin green). These plots are marginalized forms of figure 6. The distribution means are 
significantly shifted in each case, and the tails even more so. 



We now argue that there are two reasonable ways to get around this problem. One is to embrace the 
constraints, but to make them correspond to some regime in theory space that we wish to model. This 
is the approach that we are working towards when developing the parameterization introduced in [18]. 
The alternative is to evade the constraints as far as is possible by making the functions completely 
free in both k and 77. This numerically challenging approach is the one taken in [39, 40], where the 
data itself is given the freedom to choose the parameterization using a Principal Component Analysis 
(PCA) method. 

The second argument we have advanced in this brief report is as follows. Different parameteri- 
zations - even at the level of where the free functions are placed analytically - can lead to radically 
different effects on the gravitational potentials <I> and Vf, and their effects on cosmological observ- 
ables. We have focused on the the ISW effect, which gives one of the most important constraints on 
parameterized approaches. In that case significant changes to the ISW plateau can be absent due 
to cancellations in these theories, even in the case of very large modifications. The details of this 
cancellation depend on which terms are included in the extension to the Poisson and slip equations. 
Given that the goal of these parameterized frameworks is to model and detect realistic deviations 
from GR, it makes sense to consider modifications to the equations that reflect the reasonable regions 
of theory space as well as possible, as advocated in [20]. The differences highlighted in this paper 
between the parameterizations underline the delicacy of this question. 

Numerical approaches with very free functions can generate the same effect by numerically finding 
the same cancelling solutions, but models with an explicit parameterization will not in general do the 
same unless they are chosen for that specific purpose. 

In the analysis we have undertaken here the differences between the parameterizations are smaller 
than the error bars on either of them. This is a temporary situation - the data will soon improve to 
the point where the differences are significant. In the longer term as they improve further they will 
clearly pick out which part of theory space is correct regardless of which parameterization is used, 
but in the intermediate phase before then the choice will matter. Furthermore, we have focused on 
the ISW in the CMB but the arguments follow through to other cosmological observables, such as 
weak lensing, redshift space distortions and any other probe of the gravitational potentials. 

It is often argued that the purpose of these methods is not to model alternative theories, but 
simply to detect any deviation from GR+ACDM in as simple a way as possible [8, 12]. But, as we 
have shown in this paper, two different parameterizations can lead to different constraints - calling 
into question the interpretation of any individual detection if not properly put into the context of the 
class of theories that are being considered. 

Furthermore, the best way to maximize the ability of a method to find a signal like this is to 
model as closely as possible its expected characteristics, as in matched filter methods in signal analysis. 

In future work we will be extending parameterization A to model additional field components and 



modified background evolution terms. We will also be applying it to other data sets like lensing and 
cross-correlations wliich react differently to metric potentials. An inclusive approach which carefully 
includes as many phenomena as possible in a model offers the best chance of detecting any small 
deviations from GR in the upcoming era of high-precision cosmic structure data. 
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